Dynamical mean-field theory of the small polaron 



S. Chichi*, F. de Pasquale^, S. Fratini^ 1 *, D. Feinberg§ 
t Dipartimento di Fisica, Universita de L'Aquila, via Vetoio, 1-67100 Coppito-L'Aquila, Italy 
" Dipartimento di Fisica, Universita di Roma "La Sapienza", piazzale Aldo Moro 5, 1-00185 Roma, Italy 
Laboratoire d' Etudes des Proprietes Electroniques des Solides** , Centre National de la Recherche Scientifique, BP 166 3804-2 

Grenoble Cedex 9, Grenoble, France 
(February 1, 2008) 

A dynamical mean-field theory of the small polaron problem is presented, which becomes exact in 
the limit of infinite dimensions. The ground state properties and the one-electron spectral function 
are obtained for a single electron interacting with Einstein phonons by a mapping of the lattice 
problem onto a polaronic impurity model. The one-electron propagator of the impurity model is 
calculated through a continued fraction expansion (CFE), both at zero and finite temperature, for 
any electron-phonon coupling and phonon energy. In contrast to the ground state properties such as 
the effective polaron mass, which show a continuous behaviour as the coupling is increased, spectral 
properties exhibit a sharp qualitative change at low enough phonon frequency: beyond a critical 
coupling, one energy gap and then more and more open in the density of states at low energy, 
while the high energy part of the spectrum is broad and can be qualitatively explained by a strong 
coupling adiabatic approximation. As a consequence narrow and coherent low-energy subbands 
coexist with an incoherent featureless structure at high energy. The subbands denote the formation 
of quasiparticle polaron states. Also, divergencies of the self-energy may occur in the gaps. At finite 
temperature such effect triggers an important damping and broadening of the polaron subbands. 
On the other hand, in the large phonon frequency regime such a separation of energy scales does 
not exist and the spectrum has always a multipeaked structure. 

PACS numbers: 71.38.+i - Polarons and electron-phonon interaction, 
63.20.Kr Phonon-electron interactions 



I. INTRODUCTION: THE SINGLE POLARON PROBLEM 



The polaron problem is an old but not fully solved problem of solid state physics. The small polaron theory which will 
be considered here assumes a short-range electron-phonon interaction and explicitly includes the lattice periodicitytrH. 
It is therefore aimed to study systems where screening is effective. This addresses for instance the situation of a 
metal consisting of different bands, one of them being narrow enough to allow for a strong coupling to phonons. In 
fact, if the crystal can be considered as made of strongly deformable molecular-like units .with narrow-band electrons 
hopping from one to another, then the conditions for a strong polaron effect can be realized!! In realistic structures, for 
example transition metal oxides or organic metals, such units exist which provide local (oscillation) phonon modes and 
are indeed strongly coupled to well-defined electronic orbitals. Recently, much interest has been .revived into polaron 
theory, due to important classes, of materials, including the high-temperature superconductorsa and the "colossal" 
magnetoresistance manganitesQ a. In fact, in the insulating parent phase of superconducting cuprates, polarons have 
been unambiguously detected by optical measuremeiitsfl13, and some evidence of strong electron-phonon coupling 
effects has been given recently in the metallic phase&EEl. In the manganitcs, on the other hand, strong static or 
dynamic Jahn- Teller distorsions appear at the metal-insulator transitionsoS. 

The problem of a single polaron becomes relevant for low carrier density, but also by itself, as a paradigm to study 
the effect of strong coupling electron-phonon interactions. In-the intermediate and strong coupling regimes, the small 
polaron problem is already a non trivial many-body problemtlc the difficulty consists in describing the dressing of the 
electron by a coherent multiphonon cloud, moving coherently with it such as to form a quasiparticle. Pertuehatiiza 
techniques, starting either from the free electron limit or from the atomic limit (strong-coupling expansions )0'E3~t3 
fail in describing the dressing effect, ituthe intermediate regimes. On the other hand, nonperturbative solutions based 
for instance on variational ansatzEJjlLS are expected to give reliable results only for the ground state properties like 
the effective polaron mass. But no satisfactory description of the full spectral properties has been obtained so far. 

Holstein's Molecular Crystal Modelo involves tight-binding electrons coupled to dispersionless optical phonons. As 
a function of phonon frequency and electron-phonon coupling, it displays a variety of interesting regimes. The strong 
coupling regitae leads to the formation of small polarons, with a dramatic increase of the effective mass for low phonon 
frequenciesBH3. On the other hand, for exactly zero phonon frequency, an adiabatic solution can be obtained .•with 
self-trapping of the polarons appearing only above some critical coupling value, in dimensions greater than ondlJ'EjLJ. 



1 



In this context, the recent discovery of a nonperturhatiiie theory for interacting quantum problems, based on the 
limit of infinite lattice coordination (or dimensionality )EdcJ, opened a new way of attacking strong electron-phonon 
problems. A few recent papers have addressed the problem of superconductive and charge-density .wave instabilities 
of the metallic state, Q^Qse to half filling, either in weak couplingcJ, using self-consistent techniques ^30, or from the 
local impurity methocH, also including a local electron repulsion (the Holstebx-Hubbard problem)E3. Also, a solution 
at finite density in the adiabatic limit (zero phonon frequency) was obtainedu. On the other hand— exact results for 
the spectral function of a single polaron at zero temperature were recently reported by the authors^. The aim of the 
present work is to provide a complete description of the small polaron crossover based upon the knowledge of spectral 
quantities in the whole range of parameters and also at non zero temperature. The most striking features are found 
in an intermediate coupling regime where no known approximation scheme works. 

Similarly to the Hubbard model, the infinite dimensional limit, .allows to map the lattice problem onto a self- 
consistent local impurity model, here called "polaron impurity" c3c3. It consists of a single-site electron-phonon 
problem, embedded into a quantum effective medium characterized by an effective "free" propagator which has to be 
sclf-consistently determined. This mapping preserves all the complexities of the quantum dynamics of the problem, 
namely the interplay between electron and lattice fluctuations at the local level. The crucial point here is that for a 
single electron the impurity model can be analytically solved by a recursion formula for any non-interacting impurity 
propagator, leading to a continued-fraction expansion (CFE) solution for the fully dressed propagator. This unique 
feature allows to obtain at the end an exact solution for the lattice problem in the limit of infinite dimensions. This 
solution provides directly, with modest computational efforts, the ground state as well as the spectral properties in the 
thermodynamic limit and at any temperature. Thereface.it is somehow complementary to numerical sypjks performed 
in finite dimensions, such as Monte Carlo simulationsEJ'LJ or exact diagonalization of finite clustersoEj. Indeed, the 
former are limited to finite temperatures, and the latter have to deal with finite size effects. 

The main result of the present self-consistent impurity approach is that, in the crossover regime, low and high 
energy scales can be accurately described, like in the Mott-Hubbard transition problemE-X Polaron states in the low 
energy range appear as coherent strongly renormalized quasiparticle states, while at higher energies the electron is 
incoherently scattered by a quasiclassical random distortion. These features are clearly displayed in the low phonon 
frequency regime, where the spectral density displays low energy peaks coexisting with a broad and incoherent high 
energy continuum. This provides a novel and physically transparent representation of the polaron crossover. As 
the coupling strength is increased, it proceeds through successive opening of gaps in the spectral density, separating 
polaron subbands. At intermediate couplings, and decreasing the phonon frequency, one finds that the dressing of 
electron states by a multiphonon coherent cloud drastically reduces the effective electronic energy scales and leads to 
an adiabaticity "catastrophe" in the low energy spectrum. In terms of a perturbative expansion, this can be clearly 
ascribed to high order vertex dressing. It is important to emphasize that even in an intermediate regime of couplings a 
well defined polaron quasiparticle excitation is present at low energy. Also, an important feature is the appearance of 
a discrete set of frequencies where the self-energy diverges within the low-energy gaps of the spectrum. Those points 
are sensitive to disordeiH. In particular, at finite temperatures, they enhance the damping and thermal broadening 
of polaron states in their vicinity, leading to a loss of coherence of the main polaron subband. 

The paper is organized as follows . In Section II, we introduce the Holstein molecular crystal model and we discuss 
the main limiting cases in a finite dimensional lattice. In Section III, we introduce the impurity analogy and the exact 
CFE solution of the impurity model, at zero and at finite temperature. In the same section the limiting results of the 
CFE are presented and compared to the finite dimensional case. Section IV presents the general results of the CFE 
solution of the impurity problem. Section V is devoted to conclusions and to the discussion of possible extensions of 
this method to finite density, in relation with a Coherent Potential Approximation (CPA) formulation of the CFE. 

II. THE HOLSTEIN MODEL 

In this section we shall summarize the main results concerning the Holstein model in finite dimensions, in the limit 
in which analytical calculations can be performed. We especially focus on the role of dimcnsionnality, which will 
allow to discuss later how our theory (exact in infinite dimensions) can be compared with results in finite dimensional 
lattices. 

The Holstein model consists of tight-binding conduction electrons interacting with local dispersionless phonon 
modes. The Hamiltonian is 

H = - *i,j( c I,<7 C J^ + h - c -) - 9^24,<r c iA a >i + a l)+ uo^ajat (1) 

<ij>,a z,cr i 

where c\ a (ci. (J ) creates (destroys) an electron with spin a at site i, and a\(ai) creates (destroys) a phonon at site i. 
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The hopping matrix elements t^j connect nearest neighboring sites of a lattice in d dimensions and-we assume they 
give rise to a band of half-bandwidth t. This model possesses two independent control parameterstij. The first one 
is the bare coupling constant A = g 2 /uj$t = \e p \/t, where e p = —g 2 /u>o is the polaron energy obtained in the atomic 
limit (t — 0). The second one is the adiabatic parameter 7 = u>o/t. A third parameter can be conveniently introduced 
as a combination of the above ones, as a — g/uj®, with a 2 — A/7. While A and 7 are commonly used as parameters 
in the perturbative analysis, the parameter a which measures the strength of the lattice deformation involved in the 
polaron effect will show to be crucial in the strong coupling regime. Let us stress that these parameters are defined 
from the bare energy scales loq, g,t in the Hamiltonian, contrarily to usual definitions in the theory of electron-phonon 
interaction in metals where in particular the dressed phonon frequency is used. It is worth defining the following 
regimes and limits which are relevant to the Holstein model 

i) weak (strong) coupling A < 1 (> 1) 

ii) small (large) phonon frequency 7 < 1 (> 1) 
hi) multiphonon regime a 2 > 1 

iv) adiabatic limit luq = 0, finite A. 

Fig. 1 shows the corresponding regions in the (A, 7) plane. 

Throughout this paper we shall concentrate on the problem of one electron in interaction with phonons, i.e. a 
system which density is zero in the thermodynamic limit. To analyze the perturbative behaviour of the model we 
shall first discuss the simplifications due to this limit. The discussion is restricted to the zero temperature limit but 
can be easily generalized to finite temperatures (see section II B). For a single electron the Green's function in the 
site representation can be defined as 

G itj (t) = -i(0\T Cj (t)4(0)\0) (2) 

where |0) is the vacuum for phonons and electrons and the unessential spin indices are omitted. One .observes that 
there is only one possible ordering (t > 0) of the T-product, so that the function is purely retardedE3. Then the. 
standard perturbation theory is introduced in the site representation by defining the electron self-energy E^.j (o;)E3 
through the Dyson equation 

Gi,j = [Go]i,j + '^^[Go]i,k'Ek,iGij (3) 

k,l 

where Go is the free electron propagator. 

At zero density, the following simplifications hold: 

- A general self-energy diagram consists of a single electron line first emitting and then absorbing phonons. 

- The emission (absorption) of a phonon consists in subtracting (adding) a quantum of phonon frequency luq to 
the energy of the propagating electron line. 

The first statement comes from the absence of density fluctuations in the zero density limit (no bubble diagrams, no 
phonon rcnormalization). Moreover in the zero temperature limit all the phonons must be created from vacuum before 
being absorbed. To illustrate the second statement, we first notice that in a generic self-energy diagram involving 
N phonon lines, it is always possible to choose an integration contour which avoids all the poles and cuts of the 
electron Green's function, since the retarded electron propagator is analytic in the upper half plane. Then the only 
contribution comes from the poles associated to the phonon lines. 



A. Weak coupling and adiabatic limit 

The perturbation expansion of the self-energy to second order in g gives a local (k-independent) self-energy 

E 2 (w) = g 2 Go(u - loq) (4) 

where Go is the local free propagator obtained by the knowledge of tic free particle DOS N(e) as Go — J deN(e)(u> — 
e) -1 . Notice that dimensionality enters only through the free DOSE2I. The electron effective mass, in the case of a 
local self-energy, is easily calculated via 

m* dReY>(ui) , . 

= 1 Bo ( 5 ) 

m duj 
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where Eq is the ground state energy. 

Let us first consider the low phonon frequency regime. In this case for electron states lying at the bottom of a band, 
dimensionality effects enter through the band shape near the band bottom and control both the behaviour of the 
effective mass and the spectral properties.,-We assume that near the band bottom (e = —t) N(e) ~ (1 + e/i) d / 2_1 /t, 
then one has from eqs. (|],|]) for < d < 4e!3 

+ (6) 

777, 

where kd is a numerical constant.0 It is not surprising that we do not recover the expected "Migdal" result 1 + A. 
In fact this last result is obtained assuming an infinite flat band which could be the case of a metal which Fermi 
energy lies far from singularities (like Van Hove ones) in the DOS. From eq. [| we can define an effective coupling 
A = XN(E a + wo)t as in the case of Van-Hove singularities!! 2 ] in order to write m* /m = 1 + A. The effective coupling 
strength A tends to zero for vanishing phonon frequency in dimensions d > 2 (keeping A constant), while it goes to 
infinity for d < 2. Surprisingly, a perturbative analysis provides non trivial information about the adiabatic limit: for 
d > 2 we expect a free electron behaviour while for d < 2 the perturbation expansion around a delocalized solution 
fails in the adiabatic limit for any finite A. This is consistent with the non-perturbative findings of ref.E3, where it is 
shown that renormalization effecta_are absent up to a finite value of A = A c for d > 2, while in d = 1 the behaviour is 
polaronic for any finite value of Aej. 

In the opposite case of large phonon frequency, calculating the self-energy in eq. (^]) for large ujq and taking 
advantage of the asymptotical behaviour of the free propagator at high energy, it is easy to get 
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It appears as a first order expansion in power of a 2 . As_we shall see this expansion is actually resummed, in the large 
phonon frequency limit, by the use of the Lang-Firsovl!3 transformation. 

It is worth comparing the spectral properties obtained at finite bandssiidth and zero density with those derived up 
to second order in g in the classical work of Engelsberg and SchrieffeiO in the case of an infinite bandwidth. Since 
the self-energy is local one can define the spectral function as a function of energy e and frequency o£B 

A(e, u) = --Im (8) 

where e runs over the non-interacting band energies £k of a translationally invariant lattice. The electron spectral 
density Af(u>) = —(l/ir)Im,G(ui) is derived from eq. (||) by integrating over the energy distribution. 

In agreement with ref.Lj, one easily obtains a quasiparticle excitation spectrum (with ImE = 0) at energies Eq < 
w < Eq + loq and an incoherent broad spectrum at larger energies. 

In the low phonon frequency case the spectral density of the low-energy quasiparticle states can be determined by 
the low-energy properties of the free DOS. A band of coherent excitations could separate from an incoherent states, 
depending on the dimensionality. The equation which determines the band edges is 

oj - ReY,(u)) = ±t (9) 

The minus sign determines the band bottom including the ground state energy while the plus sign determines the 
band top. Using the self-energy of eq. (^) it is easy to see that near E + ujq , ReYLdi verges in d = 1 and d = 2 
as uj^ 1 / 2 and log (a;) respectively, while it is well behaved for d > 2. Uneven with ref.E3 the finite bandwidth effects 
taken into account by eq. (||) generate a gap in d = 1,2. More precisely the amplitude of the gap (in units of the 
bandwidth) between coherent and incoherent states scales with (A7) in d = 1 and with exp(— I/A7) in d = 2. In 
d = 3 the real part of the self-energy does not diverge at Eq + uio so that a sufficiently large coupling is necessary to 
fulfill eq. (Q). In this case we expect the appearance of a gap only for A greater than a certain value which depends 
on the phonon frequency and explicitly on the whole band shape. 

The discussion of this subsection suggests that, in dimensions larger than 2, perturbative expansions fail beyond 
some critical coupling above which gaps open up in the one-electron density of states. This property will be revealed 
in detail by the self-consistent local impurity theory analyzed in sections III, IV. 



B. Atomic and large phonon frequency limits 

The atomic limit is defined as the zero hopping case (t — 0). It can be understood also as an infinite coupling limit 
A —>■ 00. One considers a single electron on a single site lattice (atom) whose Hamiltonian is given by eq. (0) with 
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t = 0. In the case of zero bandwidth the Hamiltonian of eq. ([[]) can be diagonalized by the unitary Lang-Firsov (LF) 
transformation!!^ 



U = exp[aJc{a - a f )] (10) 

The effect of this transformation is to shift the phonon operators by a quantity a such as the electron-phonon 
interaction is eliminated. It introduces a new fermion, the polaron, which carries this phonon field shift 

c^Xc (11) 

where X — expa(a — a'). Once the transformation is performed, the Hamiltonian becomes diagonal and the ground 
state energy is the polaronic energy e p = — g 2 /cL>o, the excited polaron states having an energy e p + nu>o- 

Due to the presence of an electron at a given site, the lattice is deformed. The magnitude of this effect is measured 
by the local part of the static electron-displacement correlation function defined as 

Co = (m(ai + at)) (12) 

In the atomic limit one gets Co = 2a which means that the atomic ground state is that of a localized polaron i.e. 
an electron surrounded by a "cloud" represented by a coherent (Glauber) phonon state, with an average number 
< a^a >= a 2 of phonons. The electron propagator can also be calculated after the LF transformatioi£3 

GH = V^ (13) 

^ n\ uj- nu -e p 

n— 1 

The resulting spectral density appears as a Poissonian distribution of delta peaks separated by the phonon frequency 
luq. By exploiting the Lehmann representation of the Green's function one can see that such a distribution is due 
to the projection of a localized zero phonon state onto the (n-phonon) polaron eigenstates of the Hamiltonian. This 
can be usefully understood from a gedanken X-ray or optical absorption experiment, where the wavefunction of the 
localized electronic final state (with undistorted lattice) is expanded onto the (lattice relaxed) polaron eigenstates, 
which builds the electron spectral function. From eq. ( |l3| ) we see that the Green's function has a spectral weight 
close to the ground state energy which is exponentially small in the interaction strength, while the spectral weight is 
maximum for excitations involving approximately n ~ a 2 phonons. 

Let us now consider the action of the hopping. An approximation valid for large phonon frequencies is derived from 
the LF transformation, applied to the Hamiltonian with a non zero hopping term. The hopping term, modified by 
the transformation, represents the hopping of the polaron 

ti,j c i,a C j,<? ~ ¥ tijXj Xjc\ a Cj^ (14) 

The Holstein approximatiorJl consists in averaging the polaron kinetic energy on the free phonon variables, thus 
obtaining at zero temperature an effective hopping amplitude 

t itj (0\xjx j \0)=t iij e- a2 (15) 

for i,j nearest neighbors. This approximation amounts to neglect phonon emissiea and absorption during the hopping 
process. It is believed to give correct results when luq is the largest energy scalecfl. 

In the same spirit, following Alcxandrov and Ranninger, onn-can go further and use the same approximation to 
calculate the electron propagator, also for finite electron densitytj, which gives 



A(e k ,u>) = -~Im 

IT 



n—1 q 



(16) 



where e* k — exp(— a 2 ) runs over the renormalized bandwidth obtained by replacing the free hopping parameter t by 
t* = t exp(— a 2 ). This solution shows a coherent low energy quasiparticle band describing a polaron of effective mass 



2- = e« 2 (17) 
m 



located around e p , together with an incoherent structure at higher energies. 
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Let us give a physical interpretation of this result by showing that at least for the low energy states it corresponds 
to substituting the exact self-energy with the atomic one. This is valid in the case ujq >> t where a generic scattering 
process will lead electrons through intermediate states out of the band. In this scattering process the system can be 
thought as a flat band "atomic" system in interaction with high energy phonons. For frequencies near the polaron 
ground state energy e p , the atomic self-energy reads 

S(w) =Lu(l-e a2 ) + e a2 e p (18) 
Using the definition eq. (||) we get the spectral function 

A(e, u) = --Im—^ (19) 

7r uj — e* — e p 

where e* = eexp(— a 2 ) describes the renormalized band. Thus one recovers the low energy part of eq. ( fl6| ) from an 
approximation to the self-energy which is justified in the large phonon frequency regime and near the quasi-particle 
polaronic peak. However, for a 2 >> 1, i.e. when multiphonon effects are important, the validity of the Holstein 
approximation is questionable even in the case of large but finite phonon frequency. Results from small cluster 
exact diagonalizationCj show that the adiabatic ratio u>o/t must increase as a 2 to ensure the validity of the Holstein 
approximations (see fig. 1). 

Finally, let us give the result for the electron propagator in the atomic limit, at finite temperature. The Green's 
function is then defined generally by averaging on phonons only (the problem is that of a "cold" electron in a 
thcrmalized phonon bath). In the atomic limit, it is obtained in the same way as at T = 0, yielding the pole 
representatior£3 



+ oo -, 

G(u)= Y e^ 2N+1 ^ 2 I n {2a 2 [N(N + l)}^ 2 }e n ^ 2T (20) 

' lu — nu>o — e - 



n— — oo 



where N = exp{—u>o/T) is the phonon thermal weight and I n {z} are the Bessel functions of complex argument. 

Comparing this expression with eq. (|l3|), we remark that at T ^ the corresponding spectral function displays 
peaks at frequencies below the polaron ground state energy, with a spectral weight which is exponentially small at 
low temperatures. This apparent paradox of having electron states at lower energy than the ground state can be 
explained if one interprets these states as polaron states formed after absorbing n thermal phonons from the thermal 
bath, with a probability exp(— (3nu)o/2). This reduces the cost in lattice energy required to form the polaron. Since 
the polaron energy results from a balance between this (positive) cost and the (negative) electron-lattice coupling 
energy, it is possible to create states lying below the zero temperature ground state level. The price to be paid is that 
these states are incoherent, due to the incoherent (thermal) phonon distribution. Also, the chemical potential goes to 
minus infinity, allowing the fermion occupation number to be zero at any energy for one particle at finite T . 



III. THE IMPURITY ANALOGY AND THE EXACT SOLUTION FOR A SINGLE ELECTRON 

The dynamical mean field theo ry , i s developed as the exact solution of an infinite dimensionafl or infinite connec- 
tivity lattice. It has been shownEs-cll that to have a finite free electron kinetic energy the-ihopping matrix elements 
must be scaled with the square root of the lattice dimensionality or lattice coordinatior£3. A second point to deal 
with is the proper choice of the infinite coordination lattice in order to get a finite value of the ground state energy. 
In fact, a problem arises-for example in the case of an hypercubic lattice, which has a gaussian DOS with an infinite 
tail towards low energyej. For large but finite dimensions the scaling of the hopping matrix elements implies that 
the ground state energy of one electron is proportional to \fd. Therefore, the formation of a small polaron requires 
an electron-phonon coupling energy of the same order of magnitude and a coupling constant A which diverges with 
Vd. Indeed, just like the formation of a bound state from an external potential, polaron formation-by self-trapping 
requires an infinite coupling strength in an ordinarily connected lattice in infinite dimensions (see ref £3) . To overcome 
this difficulty we consider a Bethe lattice of infinite coordination d. The hopping matrix elements in_eq. ([!]) have 
been scaled as = tj1\fd with t being the half-bandwidth of the lattice. In the Bethe lattice (see ref.cZl, par. 5.3.4) 
only self-retracing paths are allowed. A restriction of the possible paths to go from one site to another allows this 
lattice to mimic a finite dimensional one also in the limit of infinite coordination, giving rise to a finite bandwidth 
semi-elliptical free DOS 

iV(e) = A V /^T^ (21) 
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which correctly simulates the low energy features of a three dimensional lattice. More generally, in interacting fermion 
problems, with this particular choice, localization phenomena can be found even in the infinite dimension limit|-4jor 
instance the Mott- Hubbard transition is correctly obtained at a finite coupling for the half-filled Hubbard modeled. 

Let us now come to the essential simplification occurring in the limit of infinite dimensions, namely the .. fa ct that 
the electron self-energy is local in space. The usual argument which holds for local electron interactionsc§E30 can be 
worked out also in the context of the Holstein model. One can carry on the standard argument but taking into-account 
also the phonon self-energies and elcctron-phonon vertices instead of the electron-electron four leg verticescj. This 

gives the scaling of the real space propagator with the intersite "Manhattan" distance R = — as G(R) oc 1/y/d . 
Using this scaling and the skeleton expansion of the vertex function^ic-can prove that the self-energy is not only local 
but that it depends only on local phonon and electron propagatori^Htj. In practice in a generic self-energy diagram 
all the internal lines are local propagators. 

Using the locality of the self-energy, the lattice propagator in the fc-space is then given by Gk(uj) = [u>— e k — E(w)] -1 
where e k is the tight-binding electronic dispersion. Writing Gu = l/^V^k^k and introducing the free DOS as 
N(e) — l/Nj^k <5( e — £ k) one has (dropping the site index for i = j) 

G(u) = [ de ^ (22) 

J w- S(cj) - e 

Notice that in infinite dimensions the properties of the lattice eater only through the free electron DOS. 

Having a local self-energy, one can demonstrate following ref.E3 the existence of an impurity model equivalent to the 
lattice problem. This can be readily seen by writing the real-space Dyson equation for the local propagator Gu (eq. 
(||)) in two steps. The first one involves self-energy contributions on sites j, with j ^= i. Once these contributions are 
resummed, one is led with a modified local propagator, noted Go- The latter can be used in the full Dyson equation 
for Gu, reintroducing the missing self-energy contributions which involve the same site i. It leads to the local Dyson 
equation for Gu — G (assuming translational symmetry) 

GH = [Go 1 (w)-S( W )]- 1 (23) 

The problem is therefore that of an impurity embedded into a medium. All the electron-phonon scattering processes 
occurring on sites other than the impurity site are contained in the effective "free" impurity propagator Go, while 
local processes at the impurity site are taken into account by the self-energy £ in eq. (|23]). 

This impurity problem can be made more physical by parametrizing it as a "polaron" Anderson impurity model 
involving a localized "d" level coupled to a local phonon, and hybridized with a fictitious conduction electron band 
" c" of dispersion E k 

H tmp = E k clc k - V k {c f k d + d f c k ) + u aU - gd f d(a + a f ) (24) 

k k 

with new impurity parameters Vk and E k being related to the propagator Go by 

G - 1 H=^- ( +0 ° dE^L (25) 

J-oo LU - E 

and 

A(£) = ^5> k 2 <5(£-£k) (26) 

k 

The original lattice Green's function it is that of the d level. Therefore, solving the problem defined by the impurity 
Hamiltonian of eqs. (|24|) for a given Go and applying the self-consistency conditions eqs. ( p2| , p3| ) one has the so called 
Local Impurity Self-consistent Approximation (LISA) which is the exact solution of a d — > oo problem. Interestingly 
enough, _Lhe above impurity Hamiltonian has been used in the past to model core level relaxation in the X-ray 
problemcl. However, in the context of the LISA approach, its significance becomes much more general. Just as the 
repulsive Anderson impurity model for the Hubbard model, it plays the role of a "paradigm" impurity model for the 
physics contained in the Holstein Hamiltonian. Though, the advantages of using an impurity parametrization of the 
d — > oo problem have been extensively reported in refEj (we also refer to the original references), we must stress, as 
a general fact, that such a parametrization is not unique in the LISA context. 
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A. The zero temperature formalism 



In the case of one single electron the Green's function at zero temperature in terms of the impurity operators, is 

G(t) = -i6(t)(0\d(t)tf (0)|0) (27) 

This function describes the propagation amplitude of the impurity electron created from the vacuum at time zero and 
destroyed at time t. Fourier transforming eq. (^) leads to the resolvent 

G{u>) = (Old 1 rft| ) (28) 

U! + 10 — H 

which has the correct prescription 5 > for convergence of the time integrals. The vacuum energy is defined here to 
be zero. 

Solving the lattice problem requires to find the solution of the impurity problem for any given Go . Let us separate 
the impurity Hamiltonian of eq. (^4|) into Hq and Hi, where Hq alone leads to the effective free propagator Go and 
Hi is the local interaction term, then an operator identity for the resolvent holds 

1111 

V^h = 7^H~ + J=T Hi 7=h (29) 

The diagonal matrix element of this operator on the impurity zero phonon state-t^|0) is the Green's function of eq. 
(p8|). To proceed further one needs to introduce the generalized matrix elemental^ 

G n , m = (0\^=d— — - — dtL^ |0) (30) 
y/n\ oj + id - H y/m\ 

so that the element Go,o will be the solution of the T — problem. 

In the case of Hi given by (pi|), to express the matrix element of the right term in eq. (|2^) one takes advantage 
of the linearity of the interaction term in the electron density operator n = Sd. Namely, introducing a set of zero 
electron p-phonon states |0,p) = (a') p /s/pl\0) one can write 

Hi = ^2d}\0,p)(0,p\d(a + a^) (31) 
p 

leading to the recursion formula for the G„ !m 's 

G n ,m — G^nSn^i Q ^ GQ n X n pGp^ m (32) 

p 

where Go n = Gq{lo — uloq) is the diagonal element of the free resolvent and X n p are the phonon displacement matrix 
elements 

X n,p = \/p+ lS ntP+ i + y/p6n, p -l. (33) 

Eq. ([32] ) is solved in matrix notation 

G" 1 = Go -1 + ffX (34) 

One immediately recognizes that, due to the particular form of X, G _1 is a tridiagonal matrix, so that the solution 
of the problem is-jjeduced to the inversion of a matrix in arbitrary dimensions. Following the lines given in ref.a (see 
alternatively ref.E3) one can express the diagonal element of the G matrix in terms of the diagonal and non-diagonal 
elements of G 1 . The local propagator (the 0, element of G) is obtained in terms of a Continued Fraction Expansion 
(CFE), as a functional of the "bare" propagator Go: 

G(u) = 2 ( 35 ) 

G °~ 1(CJ) - _ 2? 

G 1 (u - uj ) ■ ^ 

Gj\u - 2c ) - — - 9 

G (to - 3coq) - ... 
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Due to the impurity analogy, this is also the local propagator of the original lattice problem, provided that eqs. (22 2^) 
are fulfilled. As a special case, one notices that in the atomic limit, setting Gq(u>) = lo^ 1 , eq.(|35|) is nothing but 
an alternative representation of the atomic propagator, eq. (|l3|) . In the general case, the self-energy is immediately 
recognized as a functional of Go, from the self-consistency condition ( p3| ) 

E(w) = — 2 (36) 



Go 1 (lo-2lo )- 



3.9 2 



This allows to solve the impurity problem in the dynamical mean field theory. Once the self-energy is obtained for a 
given G , the local lattice propagator is calculated from eq. ( ^2|) and using eq. (p3| ) a new G is obtained. After few 
numerical iterations a fixed point Gq is reached, and the lattice local propagator G _1 = Gjp 1 — E[Gq] is determined. 
We emphasize that eq. (|3^) only involves a discrete set of values of Go and this turns out to be a drastic simplification 
of the calculation. Moreover, in the Bethe lattice the impurity and the local propagator are simply related: 

GoV)^-jGH (37) 
which, replacing eqs. (E2U23J), simplifies the calculation. 



Let us now show diagrammatically how formula (36) exa ctly sums up all of the self-energy contributions. It is 
indeed possible to relate term by term the expansion of eq. ( j36"f ) in powers of g considering Go as a parameter in the 
skeleton expansion of the self-energy. The relation is obtained through the following steps: 

i) first obtain a truncation to a given formal order g k of the skeleton expansion of the self-energy using the rules 
previously introduced. 

ii) express the internal fully interacting propagator of the skeleton expansion in terms of the "bare" impurity 
propagator Gq using the self-consistency conditions ([22 23) and expand the result to g k . The result is equal to the 



expansion of the continued fraction to the formal order g F 

It is now instructive to understand the meaning of a finite truncation of the CFE. The self-energy given in eq. (|3 
can be defined recursively 

M - 7^17 P9 \ y{p)( ; (38) 

Gq (CO — pCVo) — LW (UJ) 

where p is the stage index of the CFE for X. An A-stage truncation of the CFE is defined by neglecting E^" 1 " 1 ) in eq. 
(|38|). In the resulting diagrammatic expansion only phonon states \n) with n < N appear as intermediate states. This 
operates a selection of diagrams which is different from that based on the perturbative expansion, which by contrast 
is related to the number of interaction vertices. Indeed, at each step of the truncation an infinite set of diagrams is 
resummed, including vertex corrections (see fig. 2). As one can easily see by writing the expansion for the atomic 
limit, the parameter a 2 measures the importance of multiphonon effects and the number of phonons N needed for 
an accurate description of all the scattering processes should be much larger than a 2 . It is interesting to note that a 
self-consistent non-crossing approximation such as a Migdal scheme, always fail, since it can be put in a CFE context 



by changing all coefficients of g n in eq. ( 36 ) to 1 



B. Generalization to a thermalized lattice 

The above formalism can be easily generalized to non zero temperature. The trace performed over free phonon 
states in gives 

G(lu) = (1 - 6-*°°) e- 0n "°G nin (u) (39) 

n 

where G n ,n are the diagonal elements of the correlation matrix defined in eq. ( |30| ) and calculated by means of the 
Dyson eqs. @). 

The calculation of the diagonal elements G n<n follows the lines given in refs.CJO. The inverse of each G n ^ n {uj) is 
now the sum of an infinite continued fraction which is similar to the result at T = plus a finite fraction which 
formally takes into account the absorption processes at negative frequencies. This reads 
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where 



= G-\.)-A-B (40) 



A = ra (4i) 

(n-l)g 2 



G (uj + w ) 



1/ s (n-2)g 
Go\lj + 2u; q ) i ^ 
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.9" 



G 1 (w + nw ) 



and 



(n + l)g 2 , x 

B = —, 42 

{n + 2)g 2 V ' 



G 1 (w - w ) - 



G x (w - 2w ) - 



(n + 3) 5 2 



G 1 (w - 3w ) 



The solution of the problem now follows the same lines of the zero temperature case. 

The relation between the CFE expansion and the perturbation theory can be exploited using the simplification 
which holds in the zero density limit at non zero temperatures. The rules for constructing a self-energy diagram at 
non zero temperature are easily obtained as generalization of those obtained in section II B. 

- A general self-energy diagram consists of a single electron line emitting and absorbing phonons. 

- The emission (absorption) of a phonon consists in subtracting (adding) a quantum of phonon frequency coq to 
the energy of the propagating electron line. Associate to each process a factor 1 + /s(wo) when subtracting luq 
and a factor /b(wq) when adding wq to the electron line. 



The introduction of the temperature energy scale T/luq rules the truncation of the series of CFE eq. (p9|). For any 
finite n we also consider, in the practical evaluation of the spectral function, a finite truncation of B of size N . We 
therefore have a maximum number of phonons N + n in an intermediate virtual state of which N are emitted and 
n are absorbed from the thermal bath. Therefore the criterium of truncation valid at T = (N >> a 2 ) has to be 
supplemented by the condition n >> T/ujq. 

C. Limiting cases in the LISA approach 

We discuss here some limiting cases based upon the LISA approach and show that the main properties of the 
finite dimensional polaron problem are captured by the infinite dimensional limit. We also analyze the adiabatic limit 
which, when the polaron becomes localized, involves the breaking of translational symmetry and consequently cannot 
be achieved using the CFE, which assumes this symmetry. 

By expanding to second order the LISA self-energy given by eq. ( |36| ) and substituting the free propagator for the 
self-consistent Go we obtain the perturbative relation of eq. (H) . Then from eq. (||) and using the semi-elliptical (Bethe 
lattice) DOS of eq.(pl|) one gets the effective mass 



- = 1 - 2A 7 [1 - - 1 + 7 ] (43) 

va+T) 2 -! 



■;?) 



In the low phonon frequency regime, this becomes 



TO 



in 



1 + v/2 7 A (44) 



which shows the same behaviour as in a regular three dimensional case (see eq. . In the large phonon frequency 
regime we obtain the same result as eq. . 

The fourth order term in the self-energy expansion is 

E 4 (w) = 2 ff 4 £ 2 (w - w )£ (w - 2w ) (45) 
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^From the previously stated rules, £4 is the sum of two contributions, e.g. a fourth order non-crossing diagram (fig. 
2 b))plus a vertex correction (fig. 2 c)). In the zero density limit and in infinite dimensions, the two contributions 
are exactly equal. As a consequence, any non-crossing approximation such as the self-consistent approximation of 
Engelsberg and SchciefferL^I is not valid here. In factpthis kind of approximation can be justified by the the conventional 
Migdal's argumentfil restated for Einstein phononsEH This argument requires the condition Xuo/Ep -C 1 where Ep 
is the Fermi energyEj, which is trivially invalidated in the zero density limit where Ep = 0. 

Concerning the spectral properties most of the perturbative consideration of the three dimensional case (section II) 
are valid. Due to the fact the the real part of £ to second order does not diverge, a finite gap arises in the spectral 
density for sufficiently large values of the phonon frequency (7 > 2/(1 + 2A)). 

In the atomic limit, all the self-consistent Gq X {u} — uujq) in eq. (|35| ) must be substituted by their atomic value 
cu — uloq. The resulting continued fraction can be also independently obtained by solving directly the atomic model 
through the resolvent technique described in the previous section. The advantage of this formulation, compared to 
the LF result (eq.([l3)) is to yield immediately the self-energy as a functional of the free atomic propagator. The 
exact result at finite hopping could then be understood as the result of a Coherent Potential Approximation (CPA) 
procedure. By CPA we mean the self-consistent approximation which amcowits to substituting in the CFE form of 
the atomic self-energy a self-consistent Green's function for the atomic oneO. It is worth noting that although the 
"pole" and CFE expressions of the atomic propagator are equivalent, they give different results when extended to 
finite hopping through the CPA. For example if one extracts the self-energy for the CPA procedure from eq. JTsj), 
one recovers an expression which agrees with the exact results only for large phonon frequencies. 

The Holstein approximation described in section II B can be recovered at low energy by the CFE expansion of 
the local propagator when w is the largest energy scale in the self-consistent propagator Go- In this case, all the 
Gq 1 (uj — tuvq) with n > appearing in eq. ( |35| ) can be replaced by their atomic value giving an atomic self-energy 
which yields the exponential renormalization of the effective mass as shown in section III. 

Another instructive formula can be derived from the CFE in the limit luq — > 0. In this case all the Go(w — uluq) ia 
eq. (|||) can be replaced by Go(u>) and one recognizes the continued fraction expansion of the complex error functionEd 
which can be expressed in terms of an integral 

GM = / -^e-^l*—± (46) 

The physical interpretation of this formula is that the electron moves within a field of displacements with a gaussian 
distribution. Such distribution can be understood as the "classical" limit of the quantum probability distribution of 
local lattice displacements when the phonon frequency goes to zero while keeping the elastic energy finite. Let us 
remark that as far as the hopping self-energy term in eq. ( f46j ) is neglected one recovers a gaussian DOS as predicted 
for loo — > by the atomic limit. 

The adiabatic limit is reached when both uiq and g go to zero, keeping A fixed. In this case the CFE yields the free 
electron propagator. The problem is that one needs to consider the possibility of translational symmetry breaking. 
We have developed an independent scheme presented in Appendix A which allows an exact solution at zero phonon 
frequency keeping A finite. The ground state energy is determined by minimizing the total energy with respect to the 
lattice displacement, described by a classical variable. According to the shape of the total energy curve as a function 
of the lattice displacement (see fig. 3) we find three different regimes 

i) A < X' c . The only stable minimum corresponds to an undistorted lattice (delocalized solution, strictly free 
electron) . 

ii) X' c < A < A c . The delocalized solution is still a stable minimum, but a relative minimum appears in the potential 
at non zero lattice deformation, corresponding to a metastable localized solution (small polaron). 

iii) A > A c . The stable minimum corresponds to a localized solution, 
where A^ = 0.649.. and A c = 0.843... 

Therefore in the d — > 00 Bethe lattice a first-order localization transition occurs at A c from a delocalized free 
electron to a localized polaron. Moreover, for finite values of the coupling, the localiz&eLpolaron extends over several 
lattice shells of neighbors around a given localization site, just like in finite dimensions!^. All these features are quite 
similar to those found in regular 2-d and 3-d cubic lattices. 

To summarize the discussion of this section, examination of the various limiting regimes in the special case of an 
infinite dimensional Bethe lattice shows the consistency of this limit with a three dimensional situation. This is true 
as well in the adiabatic regime as in the perturbative and large phonon frequency regimes (shaded areas in fig. 1). 
The dynamical mean field solution, exact in infinite dimensions and presented in the following section, allows us to 
complete the phase diagram and can be thought as a controlled interpolation scheme valid at least qualitatively also 
in finite dimension d > 2. 
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IV. RESULTS FROM THE DYNAMICAL MEAN-FIELD THEORY 



Let us now turn to the self-consistent solution for the infinite-dimensional lattice with a semi-elliptical density of 



states, by solving eqs. (37, |35|). We first discuss the ground state properties, as deduced from the behaviour of S(cj) 



close to the ground state energy. 



A. Ground state properties 



The knowledge of the self-energy and of the Green's function allows to access to the ground-state properties: the 
ground state energy, which is evaluated by solving equation (||); the electropJattice local correlation function defined 
m eq. @ which can be evaluated using the Hcllmann-Fcynman theoremEj from the first derivative of the ground 
state energy with respect to g; the electron kinetic energy i.e. the average over the ground state of the hopping term 
of the Hamiltonian eq. (^) which, using again the Hellmann-Feynmann theorem, is calculated as a derivative of Eq 
with respect to t; the average phonon number in the ground state which can be obtained as a derivative of the ground 
state energy with respect to ljq. 

The ground state properties are summarized in figures 4, 5, 6, 7, 8. They illustrate the above relevant quantities 
as functions of the coupling constant A for three different values of a 2 = 1,2,5. The polaron crossover is seen as 
a continuous change from weakly dressed to quasi localized electrons. The crossover almost disappears for a 2 = 1 
while it becomes sharper for large a 2 , approaching the first-order localization transition observed in the adiabatic 
limit luq — (see Appendix A). The existence of a smooth crossover rather than an abrupt transition (for finite loq) 
corroborates the general proof previously derived by Gerlach and LowenO. However, the ground state properties are 
quantitatively calculated here for the first time for any value of the parameters. 

Let us discuss these results in more details. Concerning the ground state energy depicted in fig. 4 one sees that 
it is bounded from above by the adiabatic result (ojq — 0, a 2 — > oo) and from below by the large phonon frequency 
result (luq — oo,a 2 = 0, Eq = —A). As the coupling increases, the crossover occurs for A of the order of A c , where 
A c is the critical coupling strength obtained in the adiabatic limit. The behaviour of the effective mass is shown in 
fig. 5. For large A it increases with a 2 , but-remains smaller than Holstein's prediction of eq. (|l7|), which is attained 
asymptotically only for very large couplingsti3. 

Notice also that for a 2 = 1 we do not observe any change in the curvature of the effective mass, showing that no 
appreciable crossover occurs but rather a smooth increase of the effective mass from 1 towards exp(l) ~ 2.7. On the 
other hand, one notices that for A < A c , the effective mass diminishes as a 2 increases, in agreement with the adiabatic 
(a 2 — > oo) prediction of having unrenormalized electrons for small couplings. Therefore, the mass renormalization as 
a function of a 2 behaves in opposite ways for A < A c and A > A c . For infinite a 2 (loq — 0), m* /m jumps from 1 to oo. 

The spectral properties, discussed in the next paragraph will clarify this singular behaviour which indeed reflects 
the breakdown of the perturbation theory for A > A c . 

The ground state kinetic energy is shown in fig. 6. Again a crossover is found as a change of curvature only for 
a 2 > 1 and becomes sharper as the adiabatic limit is approached. Figure 7 displays the electron-phonon correlation 
function (i.e. the local deformation of the lattice). To make a comparison with the well defined adiabatic limit, it is 
convenient to scale this quantity by the strong coupling value 2a. A sharp crossover towards large electron-lattice local 
correlations is found for large a 2 . Similarly, the number of phonons in the ground state, shown in fig. 8 attains the 
value a 2 only asymptotically for very large coupling. As a general property, one must stress that the Holstein values 
for all the above quantities are obtained assuming a local lattice deformation. The gradual behaviour we find towards 
these values is due to the finite extension of the polaron, i.e. of the electron wavefunction and lattice deformation 
over several shells of lattice neighbours. This is also true in the d — > oo limit, and in particular causes the kinetic 
energy to be non zero, even in the adiabatic limit for A > A c (see appendix A) . 

Recently, a numerical studyE-il by diagonalization on small clusters has led to the conclusion that the polaron 
crossover occurs when both conditions a 2 > 1 and A > 1 are fulfilled. According to this interpretation, to have a 
polaron one requires that A > 1 for 7 < 1 or a 2 > 1 for 7 > 1. CFE results are in qualitative agreement with this 
statement since we observe no appreciable crossover for any value of A provided a 2 < 1, and in the opposite case the 
crossover is found at around A ~ 1 and becomes sharper as a 2 is increased. A better understanding of this behaviour 
can be gained by plotting the effective mass in the whole parameter space (7, A) fig. 9 a). The isolines corresponding 
to large effective mass define the polaron region. We see that for large A and 7 the effective mass depends only on 
a 2 as predicted by the strong coupling theory. As 7 decreases, the crossover gets sharper until it becomes a real first 
order localization transition for 7 = at the adiabatic critical value JL. Finally it must be remarked that all these 
results are in qualitative agreement with Monte-Carlo simulationscd'L 2 ], which show that at finite phonon frequency 
the ground state properties in the polaron crossover are not much dependent on the dimensionality 
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B. Spectral properties at T = 



As in other strong coupling problems such as the Hubbard model, the standard mean-field or variational techniques 
(based here on the Lang-Firsov approximation followed by phonon averaging) do not allow to go beyond the low 
energy properties. Instead, the dynamical mean-field theory provides a way to explore the whole electron spectrum. 

The spectral properties are directly extracted from the knowledge of the local propagator, and reflect the structure 
of the excited states. Formulas (||J|), valid in the case of a local self-energy, obviously apply in the case of the 
dynamical mean-field theory. 

We focus here on the behaviour of the spectral density and of the self-energy. In particular, a nonzero ImT, reveals 
an incoherent scattering due to emission and absorption of phonons. More precisely, given a band in J\[{u) of reduced 
width t* we can determine if an excitation at a given frequency has a coherent or an incoherent character. To do 
this we expand the spectral function around the peak located at to* — ReT.(uj*) = e. This relation defines a pole 
w*(e) which, in the case of a translationally invariant lattice, gives the band dispersion. Then assuming around the 
pole to* an effective mass inversely proportional to the effective bandwidth and using eq. with Eq — > to* we get 
approximately to — e — ReT,{u>) ~ (t/t*)(u> — u>*) and 

It* r 

A(e,w)zi (47) 

IT t (LU — LO*y + I Z 

where T = —(t*/t)IrnE(uj*). Then the excitation is coherent if its lifetime (1/r) is much greater than the characteristic 
time of the (renormalized) hopping processes l/t* i.e. r<f*. This gives the coherence condition JmS(w') <C t (and 
not t*). If the coherence condition holds, the quasiparticle pole u>*(e) is well defined since the spectral function has a 
sharply defined peak at lj = u)*(e) with a width much less than the renormalized bandwidth. 

The scenario for polaron formation can be analyzed from figs. 10, 11, 12 where the spectral density and the 
imaginary part of the self-energy are shown in the large, intermediate and low phonon frequency regimes respectively, 
for increasing values of A. Generally speaking, in all regions of parameters, the condition to form a (quasi)particle 
of large effective mass such as a small polaron is that a narrow coherent band emerges at low frequency. We call it 
the zero-phonon polaronic band. It is worth noticing that when a polaron with a large effective mass is formed one 
always observes several bands in the spectra, and the one at lowest energy is mxjectly coherent (JmS = 0), since this 
band lies entirely below the minimum energy Eq + loq for inelastic scatterings. Though this feature is common to 
the spectra in all parameters' regions, the way in which a polaronic behaviour shows in the spectral properties is very 
different according to the value of the adiabatic ratio 7. 

Let us first discuss the large phonon frequency regime. In this case (see fig. 10) the formation of the polaron from 
the point of view of the spectral properties is a smooth crossover. For 7 > 2, the spectra always display a multipeaked 
structure (see also fig. 9 b)) in which the n — polaron band and the edges of the first excited bands are perfectly 
coherent. These structures are subbands corresponding to polaron states with n phonons excited, and can be easily 
understood by switching on the hopping term from the atomic limit. For very large-jy, all the subbands tend to have 
the same width, in agreement with the predictions from Holstein's approximation!^, see eq. (|ilf). For intermediate 
values of 7, increasing A causes the shrinking of the width of the low energy bands and decreases their spectral weight, 
while higher frequency bands become more important. The crossover occurs for A ~ 7 (a 2 ~ 1) (cf. 12-c in which 
a 2 = 2). One notices that the envelope of the peaks in the imaginary part of the self-energy tends to reproduce the 
envelope of the bands in the spectral density, shifted by coo, as can be deduced from the CFE. Thus, for a 2 > 1, both 
the weight and the damping of the subbands tend to increase with their index and then decrease following a roughly 
Poissonian envelope, with a maximum at to ~ a 2 ujQ + E p ~ 0. 

Let us now turn to the low phonon frequency regime, in which the formation of polaronic bands is qualitatively 
different and exhibits novel features. ^From figure 12 we see that a polaronic band emerges from an incoherent band 
around A ~ 1. Increasing the value of the coupling, more and more bands emerge, having very small bandwidth. 
Notice that the low energy structures are not resolved on the scale of fig. 12 c), but they can in fact be accurately 
calculated by the CFE and are shown in detail in fig. 13. As we shall see below, the relative distance between 
polaron subbands is expected to be less than loq, due to lattice displacements effects. The inverse lifetime ImT,, on 
the opposite, reproduces the pattern of the polaron subbands at energies shifted by loq. Therefore, if the bandwidth 
renormalization is strong enough, each n-th order excited band splits further into a doublet of bands separated by 
a gap, where the main one is coherent and the secondary one is incoherent. The coherent subbands turn out to 
have equivalent heights and can be interpreted as coherent quantum tunneling out of a distorted lattice site: they 
correspond to coherent polaron bands with n — 0, 1, ... excited phonons. 

A qualitative understanding of the low and high energy excitations at small 7 can be deduced from a comparison 
with the adiabatic limit results (appendix A). Fig. 14 shows the spectrum of fig. 12 c), in correspondence with the 
adiabatic potential relative to the same value of the coupling A. It is clear from figs. 14 that a low energy scale can 
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be defined as the region where the spectrum consists of separated subbands. The separation of energy scales can be 
understood by considering the effect of finite phonon frequency i.e. by considering quantum corrections to the classical 
lattice approximation. This yields a series of low energy bands which are roughly centered around the position of the 
quantized levels of the ground-state adiabatic potential V$(X) given in eqs. ( 59|6(i| ) and fig. 14. As the energy is 



increased from Eo, hese levels become more and more hybridized with the excited adiabatic continuum associated to 
the undistorted lattice starting at V^(0) = — t. In the real spectrum, one actually observes that the very narrow 
bands merge into a broad structure just around this energy level, thus defining the amplitude of the low energy region 
as \Eq\ — t. One can also evaluate from the data shown in figs. 12 c), 13, 14 a small negative deviation of about 2% 
of the first band spacings from loq, and in general the first four narrow bands in the figure are not exactly equally 
spaced. This can be explained in terms of the adiabatic potential picture by noticing that the curvature of the total 
adiabatic energy near the distorted minima is smaller than near the undistorted position. In fact, a prediction based 
upon linearization of the adiabatic potential yields a 1.65% deviation from coq for the distance of the first excited level 
from the ground state of the adiabatic potential. 

The nature of the high energy ppJjt of the spectrum is explained using a strong coupling adiabatic approximation, 
namely the u> — > limit of the CFEo (see eq. 46). This approximation describes the broad structure as an envelope of 
resonances separated by a vanishingly small ujq. The resulting spectrum is shown in fig 14 (upper panel) and fits very 
well the high energy part of the spectral density. This fit deteriorates at intermediate energies where spiky structures 
appear, due to the effect of a finite ljq. The self-energy obtained from fl46| ) indicates that all the states at intermediate 
and high energies have an incoherent character. 

In the intermediate frequency regime, the structure of the low energy bands when these are well separated is 
curious and follows some rules that can be deduced directly from the CFE expansion (see fig. 11 d)): the lowest 
energy band is always coherent (7m£ = 0) and its shape resembles the original semi-elliptical unrenormalized DOS 
even if some asymmetry is observed towards its upper edge. This is not obtained by the usual (Holstein) strong 
coupling approximation, and could be depicted by a much larger effective mass at the top band-edge. On the other 
hand, the higher order bands acquire a complexity which can be labelled by the number of "substructures" that can 
be recognized in the band shape and increases with increasing energy. Moreover, by comparison of figures 12 and 11 
we see that the transition with increasing A from a single band structure to a multi-peaked structure at low frequencies 
is much more abrupt for small than intermediate 7, and in the former case it turns out to occur around the adiabatic 
critical value of the coupling A c . 

Let us now focus in more detail on the mechanism for gap opening in the low and intermediate phonon frequency 
regimes (7 < 2). The number of well-separated subbands is shown in fig. 9 b) as a function of the parameters A 
and 7. One notices that there is a large region of the parameters' space in which the spectral density consists of 
only two structures: a single coherent polaronic band and a high energy band which is mostly incoherent. In this 
region, approximation ( |To| ) which, apart from the low-energy coherent peak, displays an incoherent spectrum made 
of separated subbands, is not correct. The extent of this region reduces as the adiabatic ratio is reduced and for 
very low phonon frequency the system undergoes a rapid crossover to a multi-peaked structure at around A c . The 
emergence of a coherent band, separated by a finite energy gap from the continuum can be seen as A increases on fig. 
15. As a critical value of A is approached, a pseudogap appears together with a large damping of the states in the 
same energy range. At higher A, the gap is formed and keeps increasing until a second gap is formed at higher energy, 
and so on. By the way, let us underline that the total spectral weight carried by the polaron bands, when they are 
well separated, is in general different from that obtained in the atomic limit. For instance, the n = band spectral 
weight is always larger than e~ a . 

An effect related to the formation of subbands, which arises in the intermediate coupling regime, is the divergency 
of the self-energy at one and subsequently more and more frequencies located in the gaps as the coupling strength is 
increased. In the presence of disorder, external excitations or at finite temperature, a non zero spectral density can 
appear around these energies. We thus expect such excited states within the gaps to have a huge damping and to 
be localized, while the ground state can keep its delocalized character. For that motivation the authorsLJ who found 
such behaviour in the context of the Holstein-Hubbard model named it dynamical localization. We must stress that 
the relevance of such a phenomenon concerning the mobility properties such as the ac conductivity cannot be tested 
using the CFE formalism. However, we draw attention to the fact that this is a very general phenomenon, which 
occurs here for a single electron as a polaronic feature i.e. is related to the multi-peaked structure of the Green's 
function. 

The occurrence of such self-energy divergencies requires that the states at the edges of two consecutive bands have 
an infinite lifetime (JmE = 0) at zero temperature. In this case, the self-energy fulfills the eq. ([)]) at the extrema 
of the energy gap, so that the real part of id — £ changes sign within the gap without crossing the values ±t. As a 
consequence a divergency in ReT, occurs at a point ujl located in the gap. In the Bethe lattice we have the condition 
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u L - 1 -G{lo l - lu q ) - £ (2) K) = 



(48) 



where T^ 2 \uji,) is the second stage expansion of the self-energy in the CFE of eq. (|3|). Equation ( [48| ) may have more 
than one solution, i.e. many such points are expected for large values of A (see fig. 13). In the atomic limit (A — > oo) 
there is an infinite number of them, located in the gaps between the peaks of the atomic spectral density. 

In the (A, 7) plane, at least one self-energy divergency exists on the right of the dashed line in fig. 9 b). We note that 
a self-energy divergency is formed only in a preexisting energy gap. Accordingly, when increasing the coupling strenght 
a gap appears first, separating a coherent polaronic peak from an incoherent excited band, and then a self-energy 
divergency occurs in the gap when the states at the bottom of the excited band become coherent. 

This phenomenon should affect the spectral properties in the presence of disordered or at finite temperatures (see 
below). In fact, whenever there is a mechanism which is able to give rise to excited states in the proximity of those 
special points (this can be due to the lorentzian tails induced by disorder) the singularity in the imaginary part of the 
self-energy becomes a finite broadening lorentzian, which implies a loss of coherence of the neighboring states. The 
consequences on the electronic spectra at finite temperatures will be analyzed in the following section. 

C. Spectral properties at finite temperatures 

In this section we present the results obtained at finite temperature using the CFE formulation of eq. (|3^) and 
following. The analysis focuses on the intermediate regipjD where no "classical" schemes of approximation are available. 

In Holstein's original treatment of polaron motioraO, a distinction is made between transition amplitudes which 
are diagonal in the phonon number, contributing to the polaron coherent motion around the ground state energy, 
and those which are non-diagonal, giving rise to hopping-like motion. While the former decrease in the presence of 
thermal disorder, the latter are thermally activated. This allows to determine a crossover temperature T ~ OAluq 
where the polaron crosses over from coherent to hopping-like motion, which is believed to hold at strong coupling or 
in the large phonon frequency regime. 

This crossover can also be observed in the one-electron Green's function. To this purpose, we show the spectra in the 
intermediate-coupling regime for increasing T (see fig. 16). The effect is twofold: first, as is already known from the 
atomic limitEZI, (see eq. p0| ), polaron peaks appear at negative energies. Secondly, scattering by thermally populated 
phonon states causes finite lifetime effects, in addition to the zero-temperature scattering processes (see Section IV 
B), as it can be seen by inspecting the imaginary part of the self-energy. Indeed, one notices that the negative- ro 
subbands are incoherent. Moreover, at low T, a small peak appears in ImT,, close to the upper boundary of the n = 
polaron band. For larger T, some spectral weight develops in this region and contributes to gradually broaden this 
band, until the gap eventually disappears. This effect is enhanced in the vicinity of a self-energy divergency, which is 
extremely sensitive to thermal disorder. This gives rise to an inhomogeneous broadening of the main polaron band. 

By evaluating the number of coherent states (with ImT, <C t) within the low energy subbands in the intermediate 
coupling regime, one can qualitatively confirm the validity of Holstein's prediction for the crossover temperature, i.e. 
T ~ 0.3 — OAojo- On the other hand, the high energy part of the spectra is slightly smoothened by the temperature. 

A different scenario holds for other values of the parameters. Results are presented in figures 17, 18, 19 where 
the low-energy part of the electronic spectra is shown for the same parameter values as in figs. 10, 11, 12, for 
T = 0.4ciJo- One notices that in the high phonon frequency regime, the gaps between polaron bands exist even at high 
temperatures, and temperature weakly affects the overall shape of the positive n-th order polaron subbands (see fig. 
17). On the opposite, for low phonon frequencies, the shape of the polaron subbands is drastically modified (see fig. 
19): the spectral weight of each of them is roughly conserved (for not too large T), but they are noticeably enlarged, 
and consequently their height diminishes. One must underline that at those temperatures, the gaps are not destroyed 
by thermal fluctuations, apart from vanishingly small spectral density tails. 

V. CONCLUSION 

We have shown how the dynamical mean-field theory can be successfully used to solve the single polaron problem 
at any temperature. The form of the propagator as a continued fraction expansion, together with the self-consistency 
condition for the non-interacting local impurity propagator, shows that this theory yields an analytic solution of the 
problem, even if some elementary numerics is required to obtain the full spectra. We have also presented in Section 
II some results in the limiting regimes in the infinite dimensional case in order to show that the use of a semi-elliptic 
free DOS gives sensible results, in agreement with the usual three dimensional solutions. This gives confidence that 
the results from the dynamical mean-field theory are quite reliable as reflecting the actual physics in dimension larger 
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than 2. On the contrary, in the one-dimensional case, one or more gaps are always present in the spectrum, i.e. 
polarons form at any couplings. 

The properties we find for the ground state are in agreement with the conventional wisdom, in particular the absence 
of a phase transition, as soon as quantum fluctuations of the lattice are taken into account. The crossover gets sharper 
as the phonon frequency is decreased for fixed A, and becomes an abrupt transition from a free electron state to a 
small polaron state at 7 = 0. This comes directly from the first order nature of the transition in the adiabatic limit, 
clearly displayed by the direct solution. 

Beyond the ground state properties, the full electronic spectrum clarifies the nature of the polaron crossover. For 
low and intermediate phonon frequencies, as the coupling is increased, several energy gaps successively open up in the 
spectrum, leaving at low energies coherent or quasi-coherent polaronic subbands. These bands can be followed up to 
very large couplings or large phonon frequencies towards the atomic limit. Therefore they are the manifestation of 
nonperturbative processes, as is clearly demonstrated by the interpretation of the CFE in terms of a diagrammatic 
expansion. Thus, even if the ground-state properties show a continuous crossover, some higher energy features show 
a qualitative change of behaviour, the phenomenon of gap opening. 

Let us make more precise the "lattice" interpretation of these polaronic subbands. Once the self-energy E(w) 
is known, and for a given non-interacting lattice DOS, an effective dispersion relation can be found for each well- 
separated subband, by finding the poles of the spectral function eq.(J8|). This corresponds to what is commonly 
understood as polaron states, i.e. bound states of an electron with a phonon cloud, a number n of phonons being 
excited in addition. For a translationally invariant finite dimensional lattice the quantum numbers underlying the 
local (integrated) spectral function are the total momentum K of the polaronlia and the excited phonon number n. 
This obviously holds only if the subbands are coherent, i.e. only for low enough n, according to the values of A and 7. 

The low phonon frequency regime exhibits in the crossover region a coexistence of extremely narrow polaronic 
subbands at low energy and a broad featureless continuum at high energy. While the low-energy features directly 
follow from the atomic or high phonon frequency limit, the continuum is directly related to the adiabatic solution. 
The behaviour near the critical point (A = A c ,7 = 0) can be understood as follows: increasing 7, the discontinuity 
in the low energy properties becomes a sharp crossover, due to very weak coherent tunneling between small polaron 
states and quasi-free states of the equivalent impurity model. Translated into the language of the lattice problem, 
it leads to the emergence of coherent "heavy" polaron quasiparticle states ("resonances"). However, the high energy 
part of the spectrum does not reveal any qualitative change in this region of parameters. 

This behaviour demonstrates how the usual concept of " adiabaticity" , as commonly employed in metals for small 
ojq, fails in the present problem: the sharp transition at A c and the occurrence of extremely narrow polaron features 
indicates on the contrary the occurrence of an " adiabaticity catastrophe" . This is due to the relevance of high order 
vertex corrections in the perturbation expansion. 

In this work the application of the LISA approach has been limited to the one-particle propagator. In fact, no 
exact procedure has been found yet to calculate analytically for instance two particle propagators. This would allow 
to compute the (dc and optical) conductivity of polarons. On the other hand, a better understanding of the polaron 
crossover could be gained by calculating the dynamical electron-lattice correlation function. 

Let us now comment about extensions of this technique to finite electron density. As mentioned before, the CFE 
expansion keeps for the self-energy the same functional structure as that of the atomic limit. Thus it is equivalent to 
a CPA approach which could be extended at finite densities. It is however easy to show (by perturbation expansion) 
that the CPA approximation fails even at the first non-vanishing order in the density. As far as spectral properties 
are concerned, the CPA failure is particularly evident in the low energy part of the spectrum, as is well known from 
equivalent approximations for the Hubbard model. Nevertheless, our main result i.e. the coexistence of low energy 
coherent and high energy incoherent structures is a picture which should be qualitatively preserved at least at low 
carriers densities. 

Finally, although at zero density it is impossible to access the correlation functions involved in the calculation of 
the optical conductivity, let us mention that our results suggest an interpretation of the infrared response of oxide 
superconductors in the insulating phase. In these materialsO, absorption spectra for low carrier densities exhibit a 
discrete set of narrow peaks at low energies plus an incoherent background at higher energies. From our point of view, 
this could be ascribed to multiphonon excitations in the intermediate phonon frequency regime. 
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VII. APPENDIX A: THE ADIABATIC SOLUTION 



In this appendix we will solve the problem of one electron moving in an infinite coordination static Bethe lattice. 
To introduce the adiabatic limit (static lattice) one notices that the free Einstein phonon Hamiltonian can be written 
as 

H ^ = E H + \ M "l X > ( 49 ) 

i 

where Pi arc the impulses, Xi are the coordinates of the ionic motion, M is the ionic mass and loq the frequency of 
each oscillator. The adiabatic limit is achieved as M — ► oo keeping k — MuJq constant. 
The electron-phonon interaction can be written as 

H ph = -g'Y j n i x i (5°) 

i 

where the coupling constant of hamiltonian eq. ([!]) is given in terms of g' by g = g' / ^/2Mloq. The polaron energy 
e p = —g 2 /ujo = —g 2 /2k is then a well defined quantity in the adiabatic limit. 
To perform an adiabatic calculation we have to do the following steps: 

i) calculate the electronic energy for a given set of ionic deformations Xi 

ii) minimize the total energy i.e. electronic energy plus lattice elastic energy with respect to the parameters AV 
An essential relation which is useful in the adiabatic limit comes from the application of the Hcllmann-Feynman 

theorem to the ground state of the static lattice. By deriving the ground state energy of the Holstein model with 
respect to the lattice deformation Xi one obtains 

X * = TFT < n i > ( 51 ) 

^From this relation it follows that in the case of a single electron, due to charge conservation, only two situations are 
possible: 

a) delocalized solution with Xi = everywhere since the total charge density per site is zero in the thermodynamic 
limit. 

b) localized solution with some finite Xi ^ around one given site. 

Therefore, for a single electron one is restricted to study electron energies for two different classes of ionic deforma- 
tions. Notice that the results quoted above are valid at any lattice dimensionality. 

In the case of the Bethe lattice, a Dyson equation for the local propagator can be written in the adiabatic limit 

( W + A-jE G m) G m = 1 ( 52 ) 
j 

then the electronic energy can be derived from the knowledge of the adiabatic electron propagator which is the solution 
of eq. ( |52| ) for a given set of deformations {Xi}. 

In the case a) (delocalized solution) the solution is trivially the free electron propagator in the Bethe lattice. The 
ground state electron energy is then E e i = —t. 

In the case b) (localized solution) the d — > oo limit together with eq. ( |o"l| ) implies that only one site is appreciably 
distorted. Calling the center-of-deformation site it is easy to see that the nearest neighbor deformation is of the 
order 1 /d, the next nearest neighbor deformation is 1 /d 2 and so on, so that the total charge can be spread on several 
shells of neighbors even in the d — ► oo limit. The main simplification of the d — ► oo limit is then that the elastic 
energy is solely determined by the 0-site deformation, for it depends on X 2 . Consequently we have two kinds of local 
propagators: one which propagates the electron from site back to site 0, which depends upon the deformation 

= — ; — — W7, ( 53 ) 
uj + g'X - l -G 

and one which propagates the electron from and to any other site, which is free 
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In the case b) a pole in the local Green's function of eq. (m) emerges out of a band. The position of such a pole 
determines the electronic ground state energy. Using eqs. (|53| , we get the following equation 



E el + g'X - ^Re{E el - yj ' E 2 el - i 2 ) = (55) 



which solution is 



E d /t = -I for x<l/2 (56) 
E el /t = -—-x for sol/2 (57) 

where x = g'Xo/t. We see that near the zero deformation state x = the two possible solutions a) and b) coincide. 
The total potential (rescaled to the hopping energy t) we need to minimize in order to get the ground state energy is 
obtained by adding to the solution a) and b) the elastic term (which is zero in the a) case). Namely 



2 



while in the b) case 



= X Tx ~l (58) 



2 

V ad=TT~ l for x<l/2 (59) 



4A 

v ad =T7 + -T+ X for ar > 1/2 (60) 
ad 4A Ax 



the result is plotted in fig. 3. There are three different regimes determined by the following critical values for the 
coupling constant A: \' c = 0.650.. and A c = 0.844.. 

i) A < X' c The delocalized solution is the stable minimum. 

ii) X' c < A < A c The delocalized solution is the stable minimum coexisting with a metastable minimum in the 
potential V^ b \ characteristic of the localized solution. 

iii) A > A c The localized solution is the stable minimum. 

In the latter case the the delocalized solution corresponds to a continuum of unrenormalized excited adiabatic 
states, (shaded area in the lower panel of fig. 14). 

By derivatives of the ground state energy we obtain the relevant properties of the adiabatic ground state as functions 
of the coupling constant. All these functions can be expressed in terms of the derivative of the ground state energy 
Eo = V a d{X = X min ) with respect to the scaled coupling parameter x, namely A = dE /dx. 

The local electron-displacement correlation function, the electron kinetic energy are determined by deriving the 
ground state energy (Section II) 

C /2a = -A (61) 



E kin /t=-E +xA (62) 

Finally, the elastic energy is calculated as a derivative with respect to ujq and scaled to that parameter so as to get a 
finite result in the adiabatic limit 

Eph/ujQ = -xA (63) 
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FIGURE CAPTIONS 



Fig. 1 A schematic plot of the regions of parameters' space (A, 7) for the Holstcin model in an infinite dimensional 
Bcthe lattice. Below the dashed line (a 2 = 1) multiphonon processes are important. On the bold horizontal axis 
the adiabatic limit holds. The point in the upper right corner is reached at the atomic t = limit. In the shaded 
area perturbation theory (small A) or Holstein's approximation (large 7) are valid. Notice that perturbation 
theory extends its validity up to the adiabatic critical value for localization A c (see text). 

Fig. 2 The CFE expansion diagrams obtained by a truncation of the CFE at the first stage and some of those obtained 
at the second stage. Diagrams a) and b),c) represent respectively the second and fourth order perturbation theory 
terms. In this case internal propagators are assumed to be free £0 • 

Fig. 3 The adiabatic potential for a) A < A^,, b) A^, < A < A c , c) A > A c (see text). 

Fig. 4 Ground state energy versus A for three different values of the parameter a 2 — 1,2,5 (triangles, asterisks, 
diamonds). Continuous line is the adiabatic limit, dashed line the strong coupling result. 

Fig. 5 Polaron effective mass in units of the bare electron mass versus A for three different values of the parameter 
a 2 = 1,2,5 (triangles, asterisks, diamonds). Arrows mark the Holstein's approximation result exp(a 2 ). 

Fig. 6 Polaron kinetic energy versus A for three different values of the parameter a 2 = 1,2,5 (triangles, asterisks, 
diamonds). Continuous line is the adiabatic limit result. 

Fig. 7 Elcctron-phonon local correlation function scaled with the strong coupling result 2a versus A for three different 
values of the parameter a 2 — 1, 2, 5 (triangles, asterisks, diamonds). Continuous line is the adiabatic limit result. 

Fig. 8 Average number of phonons in the ground state versus A for three different values of the parameter a 2 = 1, 2, 5 
(triangles, asterisks, diamonds). Arrows mark the strong coupling result a 2 . 

Fig. 9 Both ground-state and spectral properties are summarized in the A, 7 plane, a) Effective mass isolines (from 
left to right m*=l. 1,1. 2, 1.3, 1.5, 2, 5, 20). In this picture a curve with constant a 2 is a straight line starting at the 
origin (see also fig. 1). b) The number of gaps in the electron spectral density. Near the adiabatic limit all the 
lines which separate regions of equal number of gaps collapse to A c . On the right of the dashed line at least one 
point in which the self-energy diverges appears. 

Fig. 10 Spectral density (continuous line) and imaginary part of the self-energy (dashed line) in the large phonon 
frequency regime 7 = 2, for A = a)0.08, 6)0.75, c)4.0. In this and the following spectra the energies are expressed 
in units of t. 

Fig. 11 Spectral density (continuous line) and imaginary part of the self-energy (dashed line) in the intermediate 
phonon frequency regime 7 = 0.5, for A = a)0.4, 6)1.0, c)2.0. 

Fig. 12 Spectral density (continuous line) and imaginary part of the self-energy (dashed line) in the low phonon 
frequency regime 7 = 0.125, for A = a)0.7, 6)1.0, c)2.0. The delta-like divergencies of the self-energy near the 
first four tiny bands are not reported in this figure (see fig. 13). 
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Fig. 13 Spectral density and ImS with the parameters of fig. 12. The low energy part of the spectrum is shown, 
illustrating the n-th substructures with n = 1, ..,4 (resp. 0,..,3 excited phonons). For each given n > 1, the 
coherent and incoherent states are separated by a gap. The delta peaks in JmS are revealed by adding a small 
imaginary part to the frequency u). The width of each structure has been reported using different scales to 
compare the shape of each peak. The bandwidths of the four coherent peaks are respectively from left to right: 
7.2 x 1(T 7 , 9.8 x 1(T 6 , 8.33 x 1CT 5 , 5.0134 x 1(T 4 . The width of the incoherent structures are from left to right 
6.1 x 1(T 7 , 9.67 x 1(T 6 , 8.323 x 1(T 5 . 

Fig. 14 In the upper panel is shown the spectral density (continuous line) near the adiabatic limit (A = 2,-f = 0.125) 
compared with the adiabatic strong coupling result (dashed line). In the lower panel continuous line is the 
ground state adiabatic total energy as a function of the lattice displacement (see fig. 3), dashed line represents 
the lowest excited adiabatic level which is at the bottom of a continuum (shaded area). The zero-point energy 
is omitted here. See also appendix A. 

Fig. 15 Evolution of the spectral density and imaginary part of the self-energy at low energy for growing A (upper 
panel) for 7 = 0.25 and, from left to right: A = 0.7, 0.78, 0.84, 0.9. The lower panel shows w - ReH. 

Fig. 16 Evolution of the low energy spectral density and imaginary part of self-energy with temperature in the 
intermediate regime 7 = 0.25 and A — 1.0. ^,From left to right T — 0, 0.2, 0.4, 0.8. Arrows mark points in which 
the self-energy diverges at T = 0. 

Fig. 17 Spectral density (continuous line) and imaginary part of self-energy (long dashed line) for 7 = 2 at T = 0.4. 
The spectral density at T — is shown for comparison (short dashed line). Panels a) to c) refer to the same A's 
as fig. 10. 

Fig. 18 Spectral density (continuous line) and imaginary part of self-energy (long dashed line) for 7 = 0.5 at T — 0.4. 
The spectral density at T = is shown for comparison (short dashed line) . Panels a) to c) refer to the same A's 
as fig. 11. 

Fig. 19 Spectral density (continuous line) and imaginary part of self-energy (long dashed line) for 7 = 0.125 at 
T = 0.4. The spectral density at T = is shown for comparison (short dashed line). Panels a) to c) refer to the 
same A's as fig. 12. 
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fig 1 "Dynamical mean field..." by S.Ciuchi et al . 




fig 2 "Dynamical mean field..." by S.Ciuchi et al . 




fig 3 "Dynamical mean field..." by S.Ciuchi et al . 




fig 4 "Dynamical mean field..." by S.Ciuchi et al 




fig 5 "Dynamical mean field..." by S.Ciuchi et al 




fig 6 "Dynamical mean field..." by S.Ciuchi et a] 




fig 7 "Dynamical mean field..." by S.Ciuchi et al 




fig 8 "Dynamical mean field..." by S.Ciuchi et al . 




fig 10 "Dynamical mean field..." by S.Ciuchi et al . 




fig 11 "Dynamical mean field..." by S . Ciuchi et al . 
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fig 13 "Dynamical mean field..." by S.Ciuchi et al . 
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fig 16 "Dynamical mean field..." by S.Ciuchi et al . 
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fig 17 "Dynamical mean field..." by S.Ciuchi et al . 
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